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<^ '. Abstract 

ly-v , We investigate the effect of Pauli blocking on universal two- and three-body states in the medium. 

Their corresponding energies are extracted from the poles of two- and three-body in-medium scat- 
tering amplitudes. Compared to the vacuum, the binding of dimer and trimer states is reduced by 

—* • the medium effects. In two-body scattering, the well-known physics of Cooper pairs is recovered. 

bX). In the three-body sector, we find a new class of positive energy poles which can be interpreted as 

^ , Cooper triples. 
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I. INTRODUCTION 

Ultracold Fermi gases offer a unique possibility to investigate many-body phenomena in a 
controlled environment \Mi^ ■ In dilute systems of two-component fermions, the interactions 
are characterized by the S-wave scattering length. Close to a Feshbach resonance, the 
scattering length can be tuned experimentally by varying an external magnetic field. In 
particular, the crossover from the BCS limit of weakly interacting fermions to a BEC of 
bosonic dimers by tuning through a resonance has been studied in great detail |^. The 
behavior of such a system is constrained by universal relations that involve the so-called 
contact, which measures the number of pairs of fermions with different spins that have small 
separations [sl, |6| . 

More recently, ultracold gases of three-component fermions have also been investigated. 
The interest in such systems has various motivations. First, the manifestation of the Efimov 
effect [7| has been studied in systems consisting of three hyperfine states of fermionic ^Li 
atoms. A resonant enhancement of the recombination rates at certain values of the scattering 
lengths was observed in experiment |8|, |9|. These observations were analyzed theoretically 
and traced back to the appearance of an Efimov trimer close to the three-atom threshold 
10Hl3|. Subsequently, the direct association of Efimov trimers was also achieved [ij, [l5| . 



A second line of research has focused on the phase structure of such systems [16Kl9[ . In 
these theoretical studies, two components are typically paired while the third one remains 
unpaired. This mechanism can be regarded as a generalization of the BCS case. Moreover, 
the BEC-BCS crossover has also been investigated in a three component system. In Ref. J20|, 
the dynamics of such a system was analyzed on time scales long enough to see two-body 
physics but short enough to be able to neglect Efimov states or three-body collisions. For 
three-component fermions in an optical lattice, the formation of a superfiuid phase at weak 
coupling and a "trion" phase of three-fermion bound states at strong coupling has been 
predicted |2l|. 

In this work, we combine both lines of research and investigate three-body correlations 
in the medium. We investigate the effect of Pauli blocking induced by the presence of a 
Fermi sphere on universal two- and three-body states in the medium. Their corresponding 
energies are extracted from the poles of two- and three-body scattering amplitudes in the 
medium. A similar study was carried out in Ref. [22| for the case of a fermion immersed in 
a Fermi sea interacting with two heavy bosons. The Born-Oppenheimer approximation was 
used to map the system to an effective two-body problem and calculate the dependence of 
the universal spectrum of Efimov trimers on the Fermi density in that case. In Ref. |23|, the 
modification of the Efimov spectrum for three equal-mass fermions when one of the fermions 
is embedded in a Fermi sea was calculated numerically and the modification of the universal 
scaling behavior by the background density of fermionic particles was investigated. 

Here, we investigate the medium modifications for three equal-mass fermions all of which 
are embedded in a Fermi sphere. We solve the two- and three-body scattering equations for 
this system (cf. Ref. [lO|) in the medium and present a detailed study of the poles of the 
in-medium scattering amplitude. In particular, we study the emergence of positive energy 
three-body poles analog to the Cooper pairs in the two-body system. A similar analysis was 
carried out in Refs. [2J] for the in-medium scattering amplitude of a boson and a fermion. In 
these studies, the boson-fermion Cooper pairs were found to persist for vanishing attraction. 

We consider three distinguishable non-relativistic particles of equal mass with resonant 
interactions in a Fermi sea at zero temperature. The system is assumed to be dilute, i.e. 



kpR ^ 1, where R is the range of the interaction and kp the Fermi momentum. In this 
case, the two-body interactions of the particles are determined by their scattering length 
a. We assume the two-body interactions to be resonant, i.e. \a\ S> -R. Effective range 
corrections are suppressed and can be treated in perturbation theory. Because we are at 
zero temperature, all states up to kp are occupied. For kpa ^ 1, a perturbative low-density 
expansion can be derived [25], but for k^a ~ 1 an infinite class of diagrams has to be 
resummed and one has to resort to Monte Carlo simulations or additional expansions [3|, |26 



In this study, we include only the interaction of the particles with the Fermi sea via Pauli 
blocking. These effects dominate in an expansion in the inverse number of dimensions |27l.l28( 
and determine the qualitative behavior of the system. Other effects of the medium, such as 
the excitation of particles out of the Fermi sea through scattering processes, are neglected. 
Our theoretical framework is based on an effective Lagrangian for the fermion fields \&j, 
i = 0,1,2: 

j=0 V / k=0 

where the coupling constant Qk with i ^ j ^ k parametrizes the interaction of fermions i 
and j. The term proportional to /i is a contact three-body interaction of all three fermions. 
It determines the spectrum of three-body Efimov states in the vacuum. The explicit form 
of this term will not be required for our study, since the dependence on the three-body 
term can be traded for a dependence on the cutoff in leading order calculations |29|. For 
practical calculations, it is convenient to introduce auxiliary dimer fields dk and rewrite the 
Lagrangian in the form: 

2 / V72\ 2 

£ = ^ vpt ,5, + _ U, + ^ (a,44 - I (4^.^. + ^J^l4) ) + h vl/tvl/tvl/tvi/.vl/ivl/^ . 

i=0 \ / fc=0 

(2) 
The dimer field dk describes two interacting particles i and j with i ^ j ^ k. Using the 
classical equations of motion, the equivalence of equation ([1]) and ([2]) can be demonstrated. 
This framework has been widely used to describe the universal properties of few-body sys- 
tems close to the universal limit ISOl . It has also been used as the basis for studies of the 



Efimov effect in systems three-component fermions [10|, [13 . 



II. TWO-BODY SECTOR 

A. Vacuum case 

We are now in the position to investigate the effect of Pauli blocking on universal two- 
and three-body states in the medium. We start by briefly reviewing the vacuum case. More 



details can be found in Ref. J30|]. For convenience, we set ^ = ?7i = 1 from now on. The 



bare dimer propagator derived from the Lagrangian ([2]) is simply a constant, i/Ak. The 
full, interacting dimer propagator is given by dressing the bare propagator with fermion 
bubbles, see Fig. [H It represents the exact solution of the vacuum two-body problem for the 
Lagrangian ([2]). The diagrams constitute a geometric series, which can easily be summed. 



o 



FIG. 1. Bubble sum for the full interacting dimer propagator (thick line). The double lines 
correspond to the bare dimer propagator and the single lines indicate particle propagators. 



The result can be written as 



zD,(Po, P) = 7 -, ^ , (3) 



where /(Pq? P) is the loop function for the two-fermion bubble in Fig. [H In the vacuum, the 
loop function is 

d^q i i 



«/(Po,P) 



,|<A(2vrrf + go-|(f + q)^ + .ef-go-|(f-q)' 



te 



where P = |P| and the UV divergence of the loop integral has been regulated by a momentum 
cutoff A. The cutoff dependence is absorbed into the coupling constant gk, such that all 
observable quantities are independent of A. The two-body scattering amplitude is obtained 
by multiplying the full dimer propagator with the square of the dimer-fermion coupling, 
{—igk/2y. Matching to the amplitude for scattering of particles i and j in the center of 
mass at energy E = p^, 

up')= ,/^ . =-fA(p^o), (5) 

-l/afc -tp 4 

we obtain 

Afc 1 - 2akA/n ' ^ ' 

Note that g^ and A^ are not independent at this order and all observables depend on the 
combination gf,/ A^. The renormalized dimer propagator in the vacuum can thus be written 
as 

iP)fc(Po,P).afc = ^^ l/afc - V-Po + P74 - It . (7) 

9k '- -' 

The propagator has a bound state pole at Pq = —l/a\ + P^/4 if a^ > 0. The energy at the 

pole is composed of the binding energy —l/al and the kinetic energy of the dimer P^/4. 

The total mass is 2m = 2, as expected for a dimer state. For negative scattering length, the 

pole is on the unphysical sheet and represents a virtual state. 

B. Medium case 

We now move on to medium case. In the presence of a Fermi sphere, the loop integral 
changes to 

,/(p„.P)^/ j^^mi^U^ 'Q(if-q|-M ^ (3) 



where the theta functions encode the Pauh blocking. They ensure that the intermediate 
particles can not scatter into occupied states in the Fermi sea. This introduces boundary 
conditions for the loop integrals at small momenta. Different cases must be considered. A 
summary of the calculation is given in Appendix |A1 Here we focus on the results. 

For vanishing total momentum P the boundary conditions become simple. In this case, 
the argument of both theta functions is |q| — kp. Consequently, the integration over |q| 
starts at kp and ends at A. It is evident that only the infrared behavior of the integrals is 
modified by the Fermi sea. The renormalization of UV divergences is the same as in the 
vacuum. The in-medium dimer propagator then has the form 



iDk{Po,P) 



.167r 

at 



^-^L{Po,P) 



(9) 



with 



L(Po, P = 0) = 2kp + VPo + ie Hkp - v^Pq + «e) - H^f + v^Po + le) 



(10) 



The poles of the propagator are determined by solving 



1 2kp VPo + ie 



ak 71 



TT 



An (kp - VPo + ie) - In (kp + ^/Po + 



It 



(11) 



for Pq. If Po is negative, this equation can be written as 



1 2ki 



ttk 71 



-\ — a/ I Pol arctan — ; 

7T y kp J 



(12) 



where the ie has been omitted. This equation is formally similar to Eq. (3) of Ref. |22|] for 
the binding energy of a light fermion immersed in a Fermi sea interacting with two heavy 
bosons. In this case, the Born-Oppenheimer approximation may be used and the three-body 
problem reduces to an effective two-body problem. 

In the general case, the boundary conditions are more complex (cf. Appendix |A]) . Two 
cases have to be distinguished: P < 2kp and P > 2kp. The result for the general in- medium 
loop function L{Pq, P) is: 



(a) P < 2kp: 



L(Po, P) =P/2 + kp + y/^ [In (P/2 + kp-^/^)-\n (P/2 + A;^^ + v^)] 



+ 



kl- Po- ie 



In 



P 



4 - jP' 



In (P/2 + kp-^/^)+\n {P/2 + kp + y/^) 
-In 



4-ip2 + v/^ 



(13) 
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FIG. 2. (Color online) The energy E of the dinier pole at P = plotted against the inverse 
scattering length l/a^ for a Fermi momentum kp = 0.7/Iq (solid line) and kp = l/^o (dash- 
dotted line). In addition, four selected points are marked I, II, III, and IV on the solid line. For 
comparison, the vacuum pole energy is shown by the dashed line. In the inset, the dimer pole 
energy is displayed as a function of the total momentrum P for kp = 0.7/Iq and a = Iq. Curves 
are as above. 



(b) P > 2kp: 



L{Po, P) =2kp + n^/^ + v^ 



In (P/2 -kp + y/^)+ In (P/2 + kp - y/^) 



- In (P/2 - A;^ - v^) - In (P/2 + kp + y/a) 



+ 



-kl + Po + ie 



In (P/2 -kp-^)+\n (P/2 - kp + ^) 



In (P/2 + A;^ - v^) - In (P/2 + kp + ^/^) 



(14) 



with v^ = v^^o - Py4: + it . 

We now discuss the poles of the dinier propagator in the medium. Our aim is to recover 
the known two-body physics from the viewpoint of the pole structure and then use the same 
strategy to understand the three-body sector. First, we specify our units. Since there is one 
free length scale Iq in the calculations, we express all dimensionful quantities in units of lo'- 
the energy has the unit [I/Zq]; scattering lengths [Iq] and momenta [l//o]- 

We never find more than one pole on the physical sheet in the in-medium dimer prop- 
agator. The physical conditions under which this pole can disappear are discussed below. 
In Fig. [21 the energy of the pole, E, is plotted against the inverse scattering length l/a^ at 
vanishing momentum P = for kplg = 0.7 (solid line) and 1 (dash-dotted line), respectively. 
The dashed curve represents the dimer energy in the vacuum case. For positive scattering 
length, the energy of the vacuum pole is ii^ = — 1/a^. There is no vacuum pole on the physical 
sheet if the scattering length is negative. For non- vanishing Fermi momentum, a pole with 
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FIG. 3. (Color online) (a) The energy E depicted as function of the total momentum P for 
Ofc = — 3/o and kp = 0.7//o (solid line). The vertical dotted line gives E = kp. (b) Energy E 
plotted against the Fermi momentum kp for 0^=2/9 (solid), 4/o (dashed), 6/0 (dashed dotted); 
P = 0. 



positive energy appears in the negative scattering length region. In the limit l/ctfc — ?■ —00, 
this pole asymptotically approaches the values (fc^^o)^ = 0.49 and 1 for kplo = 0.7 and 1, 
respectively. In the positive scattering length region, the pole behaves like a vacuum pole 
if the scattering length is sufficiently small. However, the corresponding binding energy is 
reduced by medium effects. Additionally four selected points are marked on the solid line: 
I, II, III, and IV. To gain deeper insight into the nature of the pole in the in-medium dimer 
propagator, these points will be further investigated below. The parameters kp and a are 
kept fixed while the momentum P will be varied. 

In the inset of Fig. [21 which corresponds to point I, the dependence of the pole energy 
on the total momentum P is shown for ak = /q and kp = 0.7 /lo (solid line). The dashed 
line shows the vacuum pole energy as before. Medium and vacuum poles have a similar 
behavior as function of the total momentum: with increasing momentum P, the pole energy 
is increased and the binding is reduced. For the vaccuum pole, this is evident from Eq. ((Tj). 
In the medium, it follows from the dominant functional dependence of the in-medium dimer 
propagator on a = Pq — P^/4 (cf. Eqs. (TT^ [HI)). Moreover, medium effects are again seen 
to lower the pole energies. 

Next, we examine the positive energy poles more thoroughly. In Fig. [3] (a), the energy 
is plotted against the total momentum for a negative scattering length a^ = — 3/o, which 
corresponds to point II in Fig. [2l The energy of the pole is positive and continuously rises 
as the momentum P is increased until the energy reaches the value kj^ = 0.49 /l^, where 
the pole disappears. This positive energy pole can be associated with Cooper pairs J3l|. 
With this interpretation, their peculiar behaviour can be understood. Assume that the two 
particles are inside the Fermi sphere. If there is no interaction, the energy of the particles 
is just their kinetic energy. In the presence of attractive interactions, the energy of the two 
particles, given by the pole energy, is lowered. Consequently, the energy gain AE is the 
difference of the kinetic energy and the pole energy. Because the maximum kinetic energy 
of two particles inside the Fermi sea is ^1^/2 -|- kp/2, the maximum energy of the pole is also 
kp. When the total momentum of the two particles becomes too large, the pole disappears. 
This property is compatible with the intepretation as Cooper pairs, whose total momentum 
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FIG. 4. (Color online) The energy of the poles plotted as a function of the total momentum P for 
a = 21q [panel (a)] and a = SIq [panel (b)] and kp = O.JIq (solid line). The dashed line shows the 
vacuum poles and the dash-dotted line gives the kinetic energy P^/4. The horizontal dotted line 
in (b) gives kp and the vertical dotted line gives 2kp. 



is commonly assumed to be zero. Remember that the energy threshold already appeared 
in Fig. 121 But in this instance a different limit was considered. The energy of the pole 
approaches the threshold kp asymptotically in the limit l/ofc -^ — oo. The energy gain AE, 
hence, decreases in this limit. But the poles never disappear and Cooper pairs can always 
be formed in this region. 

We now turn to the dependence of the poles on the Fermi momentum in the positive scat- 
tering length region. In Fig. |3] (b), the pole energy is plotted against the Fermi momentum. 
The total momentum P is set to zero and the scattering lengths are a^ = 21q, 4/o and 6/o- 
As expected, the energy is negative ed, kp = and the pole corresponds to a bound state. 
With increasing Fermi momentum, the medium effects become stronger and the binding 
is reduced. For small kp the energy rises only slowly but at larger Fermi momentum the 
energy changes rapidly, crosses zero, and becomes positive. Hence, we observe a continuous 
crossover from bound states to positive energy poles as the Fermi momentum is increased. 

So far, we could associate the left and right regions in Fig. [2] to Cooper pairs (cf . II) and 
bound states (cf. I), respectively. In between lies the crossover region. We will investigate 
this region further at the two remaining points III and IV. Here the nature of the poles 
changes as a function of momentum. In Fig. H] (a) the energy of the pole is plotted against 
the total momentum (solid line) with a = 21q corresponding to point HI in Fig. |2l By 
comparision, the dashed line shows the vacuum pole and the dash-dotted line the kinetic 
energy P^/4. For vanishing momentum the energy of the pole is extremely reduced compared 
to the vaccuum but still negative. However, in the region around P = Iq the energy becomes 
bigger than the kinetic energy. Hence, this pole can not correspond to a bound state in this 
region. For larger momenta the energy again drops below the kinetic energy and the poles 
behave similiar to vacuum poles. 

We now turn to point IV in the crossover region. In Fig. H] (b) an entirely different 
behaviour can be observerd. Again the vacuum pole (dashed line) and the kinetic energy 
(dash-dotted line) are shown. Striking are the three qualitatively different regions in this 
graph. For momenta P < 2kp the pole seems to correspond to a Cooper pair: at P = 



8 



the energy of the poles is positive, in this whole region the energy is larger than the kinetic 
energy, and the pole disappears when the energy approaches kp. In a region around P ~ 2kF 
there is no pole at all. For slightly larger momenta, the pole reappears. At first the energy is 
very close to the kinetic energy, but for larger momentum it approaches the vacuum energy, 
as expected. These pole now behaves like a bound state. 

In summary, we have related the positive energy poles at P = to Cooper pairs and the 
negative energy poles to bound states. A finite momentum P leads to an increase in the 
pole energy. In the vacuum, the additional energy is simply the kinetic energy P^/4. In the 
medium, the pole energy also increases but the dependence on P is more complicated. In 
particular, the poles can vanish and change their character. As the Fermi momentum kp 
is increased, e.g., the binding energy is reduced by medium effects. We identified the two 
extremes I and II in Fig. [2] with the BCS and BEC domains, respectively. In between there 
is a crossover region. In this region the poles change their character as a function of the 
momentum P and they can not be uniquely related to one of the two cases. Equipped with 
this qualitative understanding of in-medium two-body physics, we move on to the three-body 
amplitude. 



III. THREE-BODY SECTOR 

A. Vaccuum case 

We start by briefly reviewing the physics issues of the vacuum case and then move on 
to the medium. In the three-body system with resonant interactions, there is a universal 
spectrum of three-body bound states with an accumulation point at zero energy, called 
Efimov states [7|. The spectrum is given by Efimov's universal equation 



< + ^, = (e-2-/^")'^-"* exp [A{0/so] nl , (15) 

where the angle ^ is defined by 



a 



tnni = -a\J E^b\ (16) 

So ~ 1.00624 is a transcendental number, and k^, is the binding wave number of the state 
labelled n*. The function A(^) was first calculated in Ref. 32] and satisfies A(— |7r) = 0. 
In the unitary limit of infinite scattering length, the spectrum thus becomes geometric. The 
qualitative features of this spectrum are determined by the scattering length a, but the exact 



energies depend on the three-body interaction in Eq. (|2]) which fixes the value of k* [33[ . The 
spectrum exhibits a discrete scaling symmetry which is evident in Eq. ( !T5|) : if the scattering 
length a and the energies Eb are rescaled by the discrete scaling factor A = exp(7r/so) and 
A~^, respectively, but k* remains fixed the spectrum is mapped onto itself. If the scattering 
length dependence of one state is known, thus all other can be obtained from the scaling 



transformation. A detailed discussion of these issues can be found in Ref. [30[. Here, we focus 
on the modification of this spectrum in the medium and on possible positive energy poles in 
the three-body amplitude similar to the two-body case discussed above. As discussed above, 
we set the three-body interaction to zero in our calculation. Thus k* is proportional to the 
momentum cutoff A. The exact proportionality factor is not required for our purpose. A 
detailed study of the Efimov spectrum and the universal scaling relations in the presence 




FIG. 5. Feynman diagrams for the ferniion-dimer scattering amplitude for zero total momentum. 
Momenta p, q, k and fermion indices i,j,r,r' are assigned as in Eq. ()17p . 



of one Fermi sphere was carried out in Ref . |23[ . We go beyond this study by considering 
three Fermi spheres and exphcitly focusing on the emergence positive energy poles in the 
three-body amplitude. Preliminary results of our study were already presented in J34| . 



B. Medium case 

The three-particle scattering amplitude in the medium can be calculated by solving an 
integral equation. In order to simplify the boundary conditions given by the Pauli blocking, 
we will constrain the total momentum of the three particles to be zero. We note that the 
Fermi sea provides a special reference frame and a non-zero momentum can not be obtained 
from a simple Galilei transformation. However, we have seen in the two-body case that a 
non-zero momentum essentially increases the pole energy. Outside of the crossover region, 
the qualitative behavior remains unchanged (cf. inset of Fig. |2]). We expect the same to be 
true in the three-body case. The Feynman diagrams for three-body scattering amplitude are 
depicted in Fig. [51 Since we are interested in the three-body singularities of the amplitude, 
it is sufficient to consider the fermion-dimer scattering amplitude where the external dimer 
propagators are amputated. Because the three particles are distinguishable, we have one 
inhomogeneous and two homogeneous contributions to the amplitude as the intermediate 
dimer can be formed in two ways [10]. The integral equation for the amplitude Aij can be 
written as 



iAij{p,k,E,Ei,Ej) = - 



9i9j 



IP 



k\-kr 



4 E-E,-E 



(p+k)^ 



+ ie 



h) 



E 



9igr 

4 



] 2 

d^q iO{q — kp) 



|p + q| 



kr 



^,=0 - -|q|<A(27r) go -^q'^ + ie E - E, - go - ^{p + q^ + ie 
X iDr{E - go, g) ■ (1 - Sir) ■ iArj{q, k, E, go, Ej) , (17) 

where the momenta and particle indices are assigned as in Fig. [5] and E is the total energy. 
After setting the energies of the incoming and outgoing particles, Ei and Ej, on shell, the 
bare coupling constants are removed by defining a renormalized amplitude: 



A^Ap,k,E) = J\Zi\\Zj\Aij{p,k,E) 



(18) 



10 



where Zi is the residue of the dimer pole i in the vacuum, 

327r 



glai 



(19) 



This renormahzed amphtude has the same poles in the three-body sector as the three-particle 
scattering amphtude. We now expand the fermion-dimer amphtude in partial waves as 

oo 

A^(p, k, E) = ^(2/ + 1) {Af^)i[p, k, E]Pi{cos Ok) , (20) 

1=0 

where cos6'fc = p ■ k/(j9/c) and Pi is a Legendre polynomial. The different partial waves 
decouple and the integral equation for the /th partial wave amplitudes is 

I —Stii /■! 
i{Al^j)i\p,k,E]=-^==== / d cos OkPiicos 9k) tij{p,k, Ok, E) 

+ zV47r^ ^j^q' f\cose,Piicose,)tUp,q,e„E) 
f^o Vlflil Jkp (2vr)^ y_i 

xDr{q,m^n)ihk,E], (21) 



where 



f (.kB F) . g(|p + k|-M(l-^..) .22^ 



and 



A.(g,i?):= 



l_lL(i?-lg^g)j 
ar IT 



(23) 



is the dimer propagator without prefactors. In the vacuum only the S-wave amplitude has 
bound state poles. This remains true in the medium and we thus focus on the poles of the 
S-wave in-medium amplitude {A^j)o[p, k, E]. The technical details of the implementation of 
the boundary conditions from the Pauli blocking are discussed in Appendix |Bl In the next 
section, we present our results for the pole structure of {A^)o[p, k, E]. 

C. Results 

In this section we discuss the poles of the the S-wave amplitude {A^)o, in general for three 
different scattering lengths. In Eq. (pT]) the three-body force dependence was traded for the 
cutoff dependence, so that the cutoff determines the three-body energy in the vacuum for 
given scattering lengths. A spectrum of two states as a function of the Fermi momentum is 
shown in Fig. O Similar to the two-body case, the binding energy of each state decreases 
with rising Fermi momentum due to medium effects. Remarkable is the difference of the 
energy loss with increasing Fermi momentum for shallow and deep states. The less bound 
state disappears through the threshold while the more deeply bound state looses only about 
5% of its binding energy as kplo is increased from to 1. This behaviour of the three-body 
spectra is generic and was always observed in our calculations. 

11 



10r 



O 

C 1- 

3 



LU 



o.U 



0.2 



0.4 



0.6 



^A 



FIG. 6. (Color online) Binding energy Eb of two states depicted in dependence of kp- at = 'o for 
A; = 0,1,2 and A = 250//o. 
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FIG. 7. (Color online) (a) Energy of three-body poles plotted against kp for a^ = — l.O^o (dashed), 
qq = — 1.25/0 (dash-dotted) and oq = — 1.5/o (solid); A = 150//o and the remaining scattering 
lengths are oi = —21q and 02 = — 2.5/o- (b) The pole energy is shown as a function of 1/ao for 
A = I5O//0 (dash-dotted), A = I6O//0 (solid) and A = 170//o (dashed); kp = 0.7//o and the 
remaining scattering lengths are ai = — /q and 02 = — 0.99/o- The horizontal line gives E = 1.5 /c^. 

In Fig. [7] (a) the energy of a generic three-body pole is plotted against the Fermi mo- 
mentum for three negative scattering lengths. As in the previous case, the binding energy 
reduces with increasing Fermi momentum. Indeed, the energy goes to zero and continuously 
rises to positive values. Hence, we have found poles with positive energy. Since the total 
momentum is zero, they can not correspond to bound states. Note the resemblance between 
this figure and Fig. E] which shows dimer poles. 

To get a better understanding of these positive energy poles, we have varied one of the 
three negative scattering lengths while keeping the other two constant, see Fig. [7] (b). The 
energy rises with decreasing 1/a, but vanishes when the value of the energy becomes l.bk\. 
For different configurations of the Fermi momenta, scattering lengths, and the cutoff, we 
have always found this threshold. The accuracy of the location of this threshold reaches to 
the third (fourth) decimal place for cutoffs of the order 100 (10) /q- In order to explain this 
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observation, we draw an analogy with the positive energy poles in the two-body case. There, 
the energy gain AE is the kinetic energy minus the energy of the pole. Hence, the energy 
of the pole can not be larger than the maximum kinetic energy. In Fig. [8] configurations of 





FIG. 8. (Color online) Configurations in the Fermi sphere for two (left) and three particles (right) 
with total momentum P = 0. 

two and three particles inside the Fermi sphere are shown for total momentum P = 0. As 
discussed in the previous section, the maximum two-body pole energy is k'^p. In the case of 
three particles the magnitude of each momentum can be kp, whereas the total momentum 
remains zero. So the maximum kinetic energy of three particles inside the Fermi sphere is 
3 X A;|./2 = l.bk'^p. We conjecture that the three-particle poles belong to a state similiar to 
a Cooper pair, but built out of three particles, which we call a "Cooper triple" . In contrast 
to Cooper pairs, these Cooper triples are fermions. If the three pair scattering lengths are 
equal, the Cooper triples are SU(3) singlets. However, for different scattering lengths, the 
SU(3) symmetry is broken. 

Cooper triples also appear if one scattering length is positive. Since the energy of the 
triples is continuous in l/oj {i = 0,1,2), the region of three negative scattering lengths 
merges into the region of one positive and two negative scattering lengths at the point 
l/oj = (the other two scattering lengths are considered constant). Therefore the pole 
energy has to remain positive in the limit l/oj — )■ 0~ to obtain Cooper triples for one 
positive scattering length. For this scenario, the Fermi momentum must be sufficiently 
large. The actual value depends on the two constant scattering lengths. Hence, Cooper 
triples also occur in this region. An analogous argument holds if two or three scattering 
lengths are positive. In all three cases, we have observed Cooper triples in our calculations. 
However, it remains to be verified that the Fermi spheres assumed in our calculation persist 
in this region. 

Next, we examine which state is energetically favorable. If Cooper pairs are built in a 
three component Fermi gas, the pairs are typically formed between two components while 
the residual component remains unpaired. Therefore, we compare the energy gain of a 
Cooper triple, 1.5 kp minus pole energy, with the energy gain of a Cooper pair, kp minus 
the pole energy. The energy gain AE of a Cooper pair and a Cooper triple are compared as 
a function of one variable scattering length in Fig. [HI The remaining parameters stay the 
same as in Fig. [7] (b). Since AE of the Cooper pair depends on the scattering length, it can 
be energetically favorable to build a different Cooper pair connected to one of the constant 
scattering lengths. To account for this, we have also plotted the energy gain of the larger 
constant scattering length. We find that the energy gain of the three-particle poles is much 
larger (note the logarithmic axis). Only near the threshold for the triple, the three-particle 
AE rapidly falls off and drops below the energy gain of both Cooper pairs. This suggests 
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FIG. 9. (Color online) IS.E is plotted against the inverse scattering length 1/ao of three-body (solid 
line) and two-body (dashed line) poles; A = I6O//0, kp = 0.7//o and the constant scattering lengths 
are oi = — /q and 02 = — 0.99/o- The horizontal dotted line shows IS.E for the constant scattering 
length 02 = — 0.99/0 • 

that Cooper triples could play an important role in three-component Fermi gases in the 
continuum. 

In principle, it should always be possible to find these positive energy poles for three 
negative scattering lengths. In contrast to the two-body case, the poles do not newly emerge 
in the medium. Primarily, the poles were bound states in the vacuum which became modified 
by the medium, see Fig. [7] (a). Thus, the Fermi momentum must be large enough to obtain 
positive energy poles. This is most easily achieved for Cooper triples emerging from rather 
shallow three-body bound states in vacuum. 



IV. CONCLUSION AND OUTLOOK 

In this paper, we have examined the infiuence of Pauli blocking on universal two- and 
three-body states. First the poles of the two-body scattering amplitude in the medium were 
regarded. We were able to recover the physics of Cooper pairs and BEC-BCS crossover 
from the pole structure of the amplitude. In particular, we found that the binding energy of 
bound states decreases with rising Fermi momentum due to medium effects. In the negative 
scattering length region, positive energy poles emerge which can be identified with Cooper 
pairs. In the crossover region, the poles show a different behaviour and their nature changes 
with the total momentum. They can not be uniquely identified as bound states or Cooper 
pairs. 

We have used the same strategy to investigate the pole structure of the three-body scat- 
tering amplitude. We found that the medium effects reduce the binding of three-body states 
compared to the vacuum. This is in agreement with the findings of Ref. [23i], where the 
modification of the Efimov spectrum for three equal-mass fermions with one of the fermions 
embedded in a Fermi sea was calculated. Moreover, we found three-body poles with posi- 
tive energy. As in the two-body sector, we observed a continuous crossover from negative 
energy poles to positive energy poles as the Fermi momentum is varied. The maximum 
energy of the poles was found to be 1.5^1^. In analogy to the connection between positive 
energy poles in the two-body sector and Cooper pairs, we have interpreted this as evidence 
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for the formation of Cooper triples composed out of three particles. These Cooper triples 
are fermions and thus can not Bose condense. The energy gain of such a triple was found 
to be larger than the energy gain of the corresponding Cooper pair over a large region of 
scattering lengths. Consequently, it appears to be energetically favorable to form a triple 
instead of a pair and an unpaired third atom in this region. 

In the case of equal pair scattering lengths, the Cooper triples are SU(3) singlets. For 
different scattering lengths, however, the SU(3) symmetry is broken (cf. Fig. [7]). If the three 
scattering lengths are large, the SU(3) breaking is small since the leading corrections to the 



SU(3) limit are proportional to the inverse scattering lengths [35 . 

How these three-body correlations affect a many-body system is an open question. It 
would be interesting to extend previous studies of the phase structure of three-component 



Fermi gases [16Hl8[ to include the triples and investigate their influence. A qualitative 
picture of the many-body structure in the SU(3) symmetric limit was given by Floerchinger 
and collaborators |36| . They argue that at small density, if the scattering length is varied 
from large negative to large positive values, the BCS and EEC phases are separated by a 
trion phase of three-body bound states. At small densities, our Cooper triples must reduce 
to the trions of Ref . (36J . A related study for three-component fermions in an optical lattice 
was carried out in Ref. |2l|. Within a Hubbard model with SU(3) symmetry, a trion phase 
of three-fermion bound states has been predicted at strong coupling and a parallel to the 
baryonic phase of QCD was drawn. 

There may also be a connection to Ref. [2J| , where Boson- Fermion (BF) interactions were 
regarded in a similiar analysis. BF pairs at positive energies were found and the ground 
state was assumed to be a Fermi Gas of BF-Cooper pairs, since the pairs are still fermions. 
The interaction of three distinguishable fermions in our case could also be regarded as the 
interaction of a Cooper pair (Boson) and a fermion of the remaining type, if the scattering 
lengths are negative (BCS region) and at least two scattering lengths are different. In this 
case, we can conjecture that the ground state of the system is a Fermi gas of Cooper triples, 
which are composites of a Cooper pair and a unpaired fermion. 

Since we have only included the Pauli blocking effects from the medium, further theo- 
retical study is required. This could for example be achieved by performing Monte Carlo 
simulations of such systems similar to the two-flavour case [3]. Such a calculation would al- 
low for more quantitative predictions of the effect. In analogy to Ref. [2^, the triples might 
also lead to a new type of superfluidity in three-component Fermi systems which could be 
observed in ultracold atomic gases. For this purpose, it would be useful to calculate the 
interactions of the triples. If their interactions are attractive, they could again form Cooper 
pairs and Bose condense. Much insight would be gained if one could calculate the energy of 
such a condensate and compare it with a BCS condensate. This would allow to determine 
under which conditions such a new type of superfluidity might occur. An experimental 
test of this scenario could be carried out with ^Li atoms where mixtures of three different 
hyperfine states with tunable interactions are already available |8|, lol, |lj, [l5 . 
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Appendix A: Medium integrals 

This section gives some details of the calculation of the loop integral for the full in-medium 
dimer-propagator. The loop integral I{Po, P) is defined as follows 

i|q|<A (27r) ^ + go - 1 (£ + q) + ie ^ _ gQ _ i (£ _ q) + ig 
After a contour integration the integral simplifies to 

d'q e(|f + q|-A:^)e(|f- q| - kp) 



^I{Po, P) = t I .,^3 '' Z U_V. ■ ■ (A2) 



As already mentioned, the theta functions are boundary conditions to the integral. We 
choose P to be aligned in z-direction and switch to spherical coordinates. The integration 
still gives 27r. But the lower q boundary depends on 6, the angle between P and q. The two 
theta functions are equivalent to the following conditions: 



f-iq):=q'-Pqx + \P'-kl>0, (A3) 

f^{q):=q^ + Pqx + \P^-kl>0, (A4) 

where x = cos 9 = Pq/(Pg). The functions f^{q) are simple parabolas, whose roots are 



/+: -^±^^{x^-l) + kj,, (A5) 



/-: +^±^^(x^-l) + A;^ (A6) 

These roots do only exist for all x G [— 1, 1] if ^ < kp. Therefore we have to distinguish the 
two cases P > 2kp and P < 2kp. 

(a) P < 2kp: 

The theta functions move the lower boundary a{x) which is 



~ + ^J^{x^-l) + kl fora;<0 
The upper boundary remains unchanged. After a rescaling: ^ = s ■ kp , q = t ■ kp, 



A = A/kp and b = (Pq — ^ + ie) /kp, the integral can be written 



as 



(27r)^ [^ Jo Jsx+y/s2(x^-l)+l 0-t^ 

+ / / ^ 7-^^^^^ • (A9) 

One can easily see that the second integral merges to the first if the substitution 
X — )■ —X is performed. 
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(b) P > 2kF. 

In this case the q integration range is 



0,f 

[0,A] 



c(x) 



and 



^+ v/^,A 



for ^/l - i < X < 1 
for 



0,-f -v/^ 



and 



Px 



+ Jcix),k 



for — 1 < a; < - 



1 



with ^/c{x) := \l k% — ^{1 — x"^). Therefore the integral becomes 



''F 4 



lI{P,,P) 



iki 



(2vr) 



+ 



+ 



-sx—yjc'(x) J.2 



-dt 



b — t ■' -sx+^/dJx) b — t 



dt dx 



1— ir /-A +2 



/i-j^Jo b-t^ 



dtdx 



1 

'72 



sx-^Jdix) J.2 



t^ r t^ 

— t Jsx+y/d{x) — t 



(AlO) 



with c'{x) = 1 — s^(l — a:^). The integrals in the first and third line are equal, similiar 
to the preceding case. 



Appendix B: Integral kernel 

In this section, we discuss the calculation of the integral 

dcos9gPi{cos9g)tikip,q,9g,E) , 



(Bl) 



with 



Uj{p,q,Og,E) :-- 



e{\p + q\-kp){l-6,,) 
E — p^ — q^ — pq cos 6g + ie 



(B2) 



which is required to derive the integral equation for the three-body amplitudes in the 
medium. This type of integral appears in the inhomgeneous as well as in the homgeneous 
part. The theta function is a boundary condition on the cosOg integration: 



^(|P + q| -kp) ^p^ + 2pqcos9q + q^ > kj, . 



(B3) 



Two cases have to be distinguished. First, if |p — g| is larger than kp, the theta function is 
always fulfilled. Consequently, the integration region is [—1, 1]. If |j5 — g| < kp, the lower 
boundary will be changed. Note that the q integration begins at kp. The lower boundary 



^3 is 



COS^n 



kp — p'^ — q^ 
2pq 



(B4) 
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The angle integration for the S-wave can now be written 

/ dx— = — / dx = [ln(c — x)] , 

J„ E — p^ — q^ — pqx + le pq Ja c — x pq 

(B5) 
with c = —{E — p^ — q^ + ie) . The result is 

\p — q\ > kp '■ a = —\ 

-l[ln(c-x)]^=lflnf ^-^'-^'+^^ + " Vlnf ^-^'-^'-^^ + " ll , 

pq pq \ \ pq / V pq J J 

(B6) 
\p — q\ < kp : a = cos 6g 



pq pq \ \ pq / v pq 

(B7) 
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